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AERODYNAMICS ON A TRANSPORT AIRCRAFT TYPE WING-BODY MODEL 

V, Schmitt 


SUMMARY 

This study carried out at ONERA is based on the DFVLR- 
F4 wing-body combination. The 1/38 model is formed by a 
9.5 aspect ratio transonic wing and an Airbus A 310 fuselage. 
This configuration has been selected as the main subject of 
a GARTEUR working group (cooperation between Germany, Nether- 
lands, Great-Bri tain and Franco), 

The aim of this paper is to survey the work done by 
ONERA, After a description of the F4 wing geometrical char- 
acteristics, main experimental results obtained in the S2MA 
wind tunnel aro discussed. Both wing- fuselage interfer- 
ences and viscous effects, which are important on the wing 
due to a high rear loading, are investigated by performing 
3D calculations. An attempt is made to find their limita- 
tions. 


1 - INTRODUCTION 


The future development of the civilian transport aircraft in 
Europe, in the light of remaining competitive in the world market, 
will require more extensive studies to be conducted pertaining to 
the reduction of operating costs of such an aircraft, A special 
importance is, therefore, given to the development of methods of 
transonic calculation that are capable of handling an entire air- 
craft configuration. In this regard, considerable efforts have 
been made in Europe 111, 12 1 and in the United States 13 I , I'* 1,151. 
However, when geometric shapes become complex, the development of 
such methods and an accurate determination of their limitations also 
raise numerous and various types of problems. 


One of the problems lies in the absence of a coherent and re- 
liable set of experimental data , corresponding to a sufficiently 
realistic aircraft configuration, permitting an evaluation of cal- 
culation methods. 


♦Numbers in the margin indicate pagination in the foreign text. 



In order to solve this problem, a working group formed under 
GARTEUR (Group for Aeronautical Research and Technology In Europe), 
which brings together representatives from research institutions and 
industries of France, Germany, Great Britain and the Netherlands, 
has decided to establish such a set of data. A model was selected 
for a wing- fuselage combination designed by DFVLR using the F4 
advanced technology transonic wing and a 1/38 scale A 310 Airbus 
fuselage | <> \ <t 

This model is being tested in Europe's largest transonic wind 
tunnels and a comparison will be made of the various experimental 
results obtained. Furthermore, a comparison of these results with 
the results derived using some ten methods of transonic calculation 
will be another interesting aspect of this cooperation. 

While awaiting the conclusions of the working group, that are 
expected to be available before the end of 1982 , the purpose of this 
report is to provide a glimpse of ONERA's contribution to these 
works. 

2 - DESIGN OF THE CONFIGURATION 


The F4 wing was designed by DFVLR as a contribution to the ZKP 
(Ziviles Komponented Programm) research program started in Germany 
in 1975 with the financial support of the German Ministry for Research 
and Technology. This program, a few details of which f.r* given in | 7 j , 
lias been conducted in close collaboration between German aeronautic 
ranufac turers and DFVLR for the purpose of providing a methodical 
contribution to transonic wing technology and minimizing the risks 
resulting from an application of this technology to civilian transport 
aircraft projects. 

The conceptual design of the F4 wing is presented in detail in 
| S | and the essential points of it are presented in this chapter. 
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2,1 - Initial Conditions 

The original intention was to design a transonic wing with a 
given conventional planar shape to be integrated with a fuselage 
provided with a passenger capacity of an Airbus A310. This planar 
shape is characterized by a tapered sweptback wing presenting a 
trailing edge break in such a manner to facilitate the integration 
of the landing gear into the wing (fig.l). The main parameters aret 
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The aerodynamic conditions to satisfy are as follows i 


cruise Mach number Mo = 0.785 

(long distance) 


coefficient of lift = 0.5 

(optimum fineness ratio) 

& cruise Mach number Mo = 0.82 

(high speed) 

coefficient of lift Cz = 0,k 


12 


Thanks to transonic technology, the relative wing thickness should 
reach 1 5°r at the root and 12 to 1 3°b on the outside. This is a consi- 
derable increase compared to the A300 which has a constant relative 
thickness of 10.5^- over the entire span and a sweepback of more than 
3° at the trailing edge. 


Owing to the large aspect ratio of the wing, its design was 
first based on a careful drawing of a basic profile presenting the 
required aerodynamic characteristics. 

2.2 - Basic Profile 


The conditions of adaptation of such a profile are based on 
the simple laws of a sweptback wing with an infinite span: 
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Factor 1.2 includes all three-dimensional affacts of span Halts 
for an aspect ratio of 10* With naan values extracted froa the 
aerodynamic wing conditions indicated in the preceding chapter, the 
following conditions are derived: 

Mach number 

coefficient of lift 
S^23'^ being equal to 25 ® • 

An iterative procedure has been adopted for the design of this 
profile for which we have looked for a "plateau" type upper surface 
pressure distribution followed by an isentropic recompression, while 
the high coefficient of lift is obtained by a high rear loading* /h. 

This leads us to the definition of the DFVLR-R4 profile | y | with a rela- 
tive thickness of 13 . 5 % and the geometric characteristics of which 
are shown in figure 2. Shown also, for illustration, (fig, 3) are 
the pressure distributions obtained in the vicinity of an adaptation 
point by calculation |10] and in the TVB wind tunnel of the DFVLR* 

2*3 - Geometric Characteristics of the Wing 

Based on this profile which aust assure the desired type of 
recompression, the wing design is formed with the following objectives: 




- 0.73 
= 0.65 


-to obtain a minimum induced drag by an elliptical lift distribution 
along the span; 

-to obtain along the wing surface a network of isobars parallel to 
the lines at X/c s constant. 
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To accomplish this, the wing tvist lav has been determined 
using a reverse method of vortex network Till with an elliptical 
lift distribution by starting with a load distribution corresponding 
to that of profile R4 at the adaptation point. The external part of 
the wing (beyond the trailing edge break) is obtained by projecting 
the basic profile in the wind bed (by the cosine of sweep angle 
23 "!) and by its adjustment to the calculated twist angle • 

For the definition of the profile at the root, iterative calcula— 
tions have been performed around the entire configuration using a 
transonic method of small disturbances I 12 I by applying the above- 
mentioned concept of straight isobar lines. 

Finally, the DFVLR-F4 wing is defined geometrically from 4 
definition sections presenting the following characteristics: 

c 'c £ ^ 

0,15 A. 5° 

0,122 1 , 8 ° 

0,122 0 , 9 ° 

0,122 - 0 . 5 ° 

The relative thickness laws and the wing twist are shown in 
figure 4, In order to obtain a simple geometry, a linear interpo- 
lation procedure has been adopted between the 4 design sections. 

This wing is integrated with the fuselage in low position with 
a dihedral angle of 4,8°, 

3 - EXPERIMENTAL STUDY IN THE S2MA WIND TUNNEL 

This study was preceded first by a 1/38 scale half-model 
( •* ;•’ 0,5877 n)_ in the GBttingen wind tunnel of DFLR (see ref, 

then by the first experimentation of the complete model (same scale) 
at NLR in the HST wind tunnel. It should be noted that the construc- 
tion of the full model and the execution of this test correspond to 
GARTEUR operations. 
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As with the case of the HST wind tunnel y the full Model was 
mounted at S2MA on a Z-shaped sting using a 6 component balance* 

The right wing is equipped with 252 pressure inlets distributed into 
7 measuring sections, plus 44 pressure inlets placed along a vertical 
plane to the fuselage. Moreover, the left wing received an instru- 
mentation for unsteady measurements (l root gage, 2 accelerometers, 

7 "Kulite" pressure transducers) and 14 colored fluid transmitters 
for buffeting investigations. 

The tests were performed at the maximum perforation rate 
of the test section (6%), Under these conditions, the wall effects 
are known to be very small. Considering its dimensions, at the 
cruising point of model F4, we can expect to have corrections of 
about : 

r + 0.003 . M . 

in Mach number 

- + 0.03° 

- in incidence 

Nevertheless, in order to meet the abovementioned quality 
conditions, corrections are made of wall effects, of the pressure 
gradient in an empty test section and of the sting effect. 

In order to assure an interference of shock wave - turbulent 
boundary layer, which is considered to be indispensable, it is 
necessary to start the transition on the upper surface and on the 
undersurface of the wing using the nominal Reynolds number of these 
tests ( <i» c . = 3.10*’, and the aerodynamic reference chord 

c = 0.1412m). The next chapter describes a few problems connected 
with this transition initiation. 

3.1 - Evaluation of the Effects Caused by the Initiation of the 
Transition on the Wing 

The transition initiator retained for all tests of model F4 
is a carborundum band 3 mm wide which has the same position on the 
upper skin and on the lower skin for all tests. On the top skin, 

it is selected as a function of the shock wave positions at the 
adaptation point and the surrounding area. It is the result of a 
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compromise between a position too far forward and a position too 
far backward which would prevent the formation of a homogeneous 
turbulent boudnary layer upstream from the shock* On the under- 
surface, this position is determined as a function of the pressure 
gradient in order to avoid too much reduction in the rear loading 
caused by a thickening of the boundary layer* The selection made 
for the F4 wing is illustrated in figure 5» 

We now have to determine the optimum size of carborundum for 
initiating the transition on the band without excessive thickening 
of the turbulent boundary layer in order to avoid a parasitic drag* 
This operation is performed in an industrial wind tunnel using vis- 
ualisations of the transition by a sublimation technique* This in- 
vestigation is usually conducted with the cruise Mach number for a 
coefficient of lift Cz in the vicinity of the appearance of bufet- 
ting. A critical analysis of the selection performed at S2MA has 
been tried at DERAT by a purely two-dimensional approach using an em- 
pirical criterion of optimum initiation The result is shown 

in figure 6 for a measuring section situated at -Zy/S = 0,636 
at the upper surface and at the undersurface in the case of cruise 
conditions. The roughness heights seem to be close to their optimum 
values. The excess thickening of the boundary layers is therefore 
evaluated to be very small* 

Nevertheless, considering the rather low Reynolds numbers that 
can now be obtained on transport aircraft models, more effort is 
needed to develop a more reliable technique of transition onset in 
a wind tunnel and to improve the prediction means in this regard of 
the device used on the created turbulent boundary layer* 

The general consequences of transition initiation may be brought 
to light by a comparison with tests of natural transitions* For 
illustration, figure 7 shows the pressure fields on the wing in 
natural transition ( f\- ) and the initiated transition ( r o ) 
for Mo = 0*75 and = 0,5, which corresponds to one of the 3 

cases of calculation selected by the GARTEUR research team which 
we will return to in chapter 4. Notwithstanding the noticeable 
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reduction of the rear loading at the external wing undersurface in /6 
TD , the shock wave position at the upper surface seems to be hardly 
affected by the initiation. On the other hand, the variation of 
the coefficient of pressure at the trailing edge of the section 
situated at - o, •'*(><) at this Mach number and also at the 

cruise Mach number (fig. ) diverges at C that are much lover in 
TC than in TN. Since this divergence is expressed by the appear- 
ance of a separation at the trailing edge, the Cz appearing with 
buffeting is also lower in TD than in TN. Likewise , we can see in 
figure 9 that the fineness ratio, at the cruise Mach number, under- 
goes a loss of about 10# at the cruise Cz relative to the configur- 
ation in TN. 

3.2.- Aerodynamic Performances 

To have a general idea of the aerodynamic quality of the con- 
figuration under study (in TD) , the drag divergence limit should 
first be studied in a plane Cz, Mo. This requires the study of the 
drag coefficient variation Cx as a function of the Mach number at Cz 
(figure 1 0 ) . In passing, it may be observed on this figure the ab- 
sence of a noticeable precritical drag. The determination of the 
drag divergence from these curves has been performed using the usual 
criterion 

^v- x /o> bf^ = 0,1 

The result obtained is shown in figure 11 and shows in particular 
that the target cruise point for the flight at an optimum fineness 
ratio is situated below this limit that physically expresses the 
presence of a supercritical flow with the appearance of shock waves. 

Finally, in order to evaluate the induced drag, it is convenient 
to examine the experimental curve 

C\ - } 
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The slope of this curve relative to the optimum theoretical value 

<Pc\ __ 

Tc£ 77 A 

is a quality index. Figure 12 presents the result obtained at the 
cruise Mach number, the slope ratio being £ - 0.882. 

The aerodynamic quality of a project may also be expressed in 
the Br4guet formula relative to the distance that the term can cov- 

er , /y 

«- ' v_ y 'o 

the maximum of which is researched for the cruise point. A plotting 
of lines at the same level, extracted from the set of stress meas- 
urements, in plane Cz, Mo (figure 13) illustrates that this condi- 
tion is quite satisfactory in the present case. 


One last point concerning the limit of a buffeting appearance is 
as follows: considering the model equipment used, several types of 
measurements allow us to make an evaluation of this limit: 

-steady measurement using the balance: c £ c fco()\ the appearance of 

buffeting is determined at = 0 , 1 ' beyond the linear variation 

of this curve, 

-unsteady measurements ( g'yS values within a given frequency L 1 

band) by the balance (dynamic roll c /i -c , a bending gage 
at the root and an accelerometer integrated in the wing. In these 
cases, which correspond to various structural responses, the ap- 
pearance of bueffeting is diagnosed by the signal divergence, i.e. 
a sudden increase in its intensity when the incidence increases 
beyond a certain value. 

Figure 14 shows that the range formed by all limits thus found 
is very wide and its interpretation a not easy. Nevertheless, the 
limit deduced from the curves of lift are likely to be pessimistic 
in a case such as this one for which the important viscous effects 
affect the flow on the airfoil. 
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If a first problem in the study of buffeting is an accurate 
determination of the appearance of this phenonmenon, an even great- 
er problem is quantification. An attempt is being made to implement 
the very simple, but also very controversial Mabey method ^[14). 

4 - PREDICTIONS BY CALCULATION AND COMPARISON WITH EXPERIMENTS 


The conditions for the calculations to be performed on the 
configuration under study vere suspended after the experiments on 
the model at NLR vere completed. The selection made on cases in 
the vicinity of the cruise point exempt from pronounced local, separ- 
ations: i 

L 2 

°.75 0.75 

Mach number Mo 

coefficient of lift Cz 

or incidence (y , 0,10 0,84 

In order to account for deformations subjected by the airfoil 
under aerodynamic loading, a minor correction of the twist lav 
(chapter 2 . 3 ) has been made as a result of a theoretical evaluation 
at the NLR. 

4.1 - Supercritical Flow on the Wing 

Calculations have been performed at 0NERA, as is the case with 
most participants, by processing the airfoil separately as an "exposed" 
airfoil. These are perfect fluid calculations using a relaxation 
method to solve the equation of full potential in a non-conservative 
form. This method isddescribed in detail in ' [ 15 ] . The code of 

V 

calculation has, nevertheless, been constantly improved and it now 
exists in a more operational form. 

Since a calculation - Cz experiment comparison is likely to under- 
go considerable effects of compensation, the calculations have been 
performed at a given incidence. This method of procedure does not 
favor a comparison with the experiment, but it does make it possible 
to evaluate important viscous effects that are expected in view of 
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the high rear loading of the wing 
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Calculations are performed using a boxed meshing technique, 
the finest mesh of which is 1 30 x 32 x 30» or 124,800 points, while 
the flow on the wing is described by 20 sections along the span with 
84 points in each (42 at the upper and lowersurface ) • The meshing, 
of parabolic type, develops spanwise. Figure 15 shows details of 
the sections closest to the root and tip of the F4 wing. It may be 
seen that the sheet outlining the slip-stream straightens as we move 
downstream from the wing. 

The variation of the maximum residue along the iterations and 
the convergence of the lift are shown in figures 16 and 17 for the 
case of the first calculation. The last figuro indicates that the 
result obtained has not completely converged, in spite of a I 85 nin 
computer time on the Cyber 170-750 computer. 

The solution obtained is shown by iso-Mach lines at the upper 
surface (fig. 18), then by the pressure field viewed in perspective 
along the upper surface and on the undersurface (fig.19)» On the 
upper surface, we may see the presence of a widely extended super- 
critical zone with a high plateau preceded and followed by relatively 
important recompressions and, on the undersurface we may see the 
high rear loading, particularly on the external part of the wing. 

The comparison of these pressure distributions with the exper- 
iment on figure 20, however, is not satisfactory. It is not satis- 
factory in the fine description, with in particular the absence of 
a definite shock wave and a too high rear loading in the calcula- 
tion, or at the level of the coefficients of local lift » 

the theoretical values of which are clearly too high. The span 
load distribution (fig.2l) confirms the high overestimation of the 
calculations, although the appearance of the curves are quite sim- 
ilar. Furthermore , the lack of agreement in the calculation-exper- 
iment comparison appears also in the other two cases treated. 

It should be remembered that the transonic calculations per- 
formed do not include the presence of the fuselage and the viscous 




j effects. In the next chapter we will try to evaluate the respective 

role of each of the two factors. 

4.2 - Evaluation of the Influence of the Fuselage and Viscous Effects 

| 

In order to determine the respective influences of the fuselage 
and of the viscous effects on the airfoil, the case of a subcritical 
flow was prepared. The following conditions were selected t 

Mach number Mo = 0,60 

incidence <X = o,7° 

In order to evaluate the influence of the fuselage, we used a 
method of singularities developed by Aerospatiale fi6] in which the 
effect of thickness is determined by a surface distribution of 
sources, whereas the lift effect is obtained using a doublets distri- 
bution on the outline of the airfoil. 

The calculations performed are, of course, concerned with the 
wing- fuselage combination, the quantification of which is shown in 
figure 22, then with the fuselage separately and the wing separately. 
Two calculations were performed for the latter: 

-one for the wing extrapolated to the symmetry plane, 

-the other for the exposed wing as adopted for transonic calculations. 

First, the wing to fuselage interferences are shown by comparing 
the pressure distributions (fig. 23) obtained in the symmetry plane 
for complete system configurations and for the fuselage separately. 

An excellent agreement between calculations and experiment may be 
noted . 


The fuselage to wing interferences are shown in figure 24, where 
a comparison is given, first, of the span loading distributions cal- 
culated on the entire configuration and on the extrapolated wing. 

It may be observed that the theoretical lift of the extrapolated 
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wing overshoots that of the full configuration by sons 8 %. 

The calculation-experiment comparison of the local lift distri- 
bution on the full system brings to mind that of the supercritical 
case (fig. 21 ). The calculated lift exceeds the experimental value 
by 30 %. 


Another concern is to bring to light the benefits of exposed 
wing calculations compared to full configuration calculations. The 
answer may be found in the results presented in figures 25 and 26 . 

By comparing the lift distributions obtained by the method of singu- 
larities, we may detect, at the wing level, relatively few differences 
between the full configuration and the exposed wing which are ex- 
pressed by a lift deficiency in the case of the exposed wing of 
about 3.5%. Moreover, the result of a calculation of the exposed 
wing using the transonic method agrees perfectly with the result 
given by the method of singularities. The pressure distributions 
calculated on the full system and on the exposed wing by using both 
methods confirms the idea that the concept of the exposed wing in 
the calculations is acceptable, at least in the subcritical case. 
Furthermore, a few published results 1 17 1, | 1H | tend to show that 
this observation is also true for supercritical flows. 

To bring to light the influence of viscous effects on this con- 
figuration, a two-dimensional approach has been tried by using 
a section at mid-span from the wing. The calculation code used 
solves the entire equation of potential in connection with the 
evaluation of the boundary layer using an integral method f 191 
In an iterative calculation process, the geometry is corrected of 
displacement effects, with the convergence intervening a few steps 
later. 


For the example treated here (fig. 27) t the two-dimensional 
calculation was performed by beginning with the perfect fluid of 
the coefficient of local lift determined by three-dimensional cal- 
culation. The agreement of the pressure distributions thus calcu- 
lated with the experiment is good, except at the undersurface at 
the rear loading level. Compared to the experiment, the calculation 
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of the perfect fluid leads to an overestimation of the local lift 
by 26$, that of the viscous fluid leads to on overestimation of 
only 4$. This proves that the influence of viscous effects on the 
F4 system is considerable and explains the poor quality of direct 
comparisons of perfect fluid calculations and experimentation. 

5 - CONCLUSION 

The DFVLR-F4 model, designed within the framework of the Ger- 
man ZKP research operation, has been selected by a research team 
of GARTEUR as a support for transonic studies of a wing- fuselage 
system used to represent a civilian transport aircraft. Experimen- 
tal results are expected to be obtained in four of the larges trans- 
onic wind tunnels of Europe to be used as a basis of evaluation of 
operational transonic calculation methods or of those being devel- 
oped. 


Ve may at this time conclude from results already collected in 
the S2MA wind tunnel that the aerodyanmic conditions of the project 
are satisfactory. However, the initiation of the boundary layer 
transition on the wing, which is imperative for obtaining turbulent 
boundary layer - shock wave interferences on the top skin using 
Reynolds number - 5 . io ( ' results in a few losses from the point 

of view of performances. 

Furthermore, the calculation - experiment comparisons performed 
at ONERA stress the presence of important viscous effects on this 
configuration. On the other hand, the influence of the fuselage 
does not seem to be a major obstacle for a comparison of experimental 
results with calculation results obtained on only the wing. 
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Fig. 1 - Planar Shape of the DFLR-F4 Wing. 


Geometric parameters: 

aspect ratio 
taper ratio 
leading edge sweep 
fuselage radius 
trailing edge break 


A =9.5 
A =0.3 
f *27.1° 
W *0.1263 
2fc/b *0.4 



Fig. 2 - Geometric Characteristics 
of the Basic Profile. 
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Fig. 18-1 so-Mach Lines Calcu- 
lated in the Vicinity 
of* a Crossing Point. 


Fig. 1 6 -Maxi mum Error Variations. 

Ft Mo=O f 75 a =0,1° 
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Fig. 18— Field Pressure Calculated in the Vicinity of 
the Crossing Point. 
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DFVLR-F4 MODEL M>=0,75 (1*0,1° 



Fig. 20-Pressure Distributions (Theoretical and Experimental) • 
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Fig. 21-Span Lift Distributions. 
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Fig, 22-Subcritical Calculations of the 
Geometric Quantification, 


Wing- Fuselage F4 Model 
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Key* a- (Singular! ties) cal- 
culation; b-Full system; 
c-Wing; d-Test; e-Full 
System. 
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Fig. 24-Effects of Fuselage-Wing Integration. 
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Key: a-Full system (singu- 
larities); b-Wing part; 
c-Exposed Wing; 
d-Singularities. 
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Fig. 25-Influence of the "exposed wing" 

design on the span lift distribution. 
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